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Abstract 

Purpose: To compare MRI sequences for breast density measurements on a 3T MRI system using IDEAL (Iterative 
Decomposition of water and fat with Echo Asymmetry and Least squares estimation) as possible physiology-like reference. 

Materials and Methods: MRI examination was performed in 48 consecutive patients (mean age 41, years; range, 35-67 
years) on a 3.0T scanner and 46 were included. All (fertile) women, were examined between days 5 and 15 of their menstrual 
cycle. MRI protocol included: T1 -turbo spin-echo (T1-tSE), T2-turbo spin-echo (T2-tSE), VIBRANT (Volume Imaging for Breast 
Assessment) before and after injection of contrast media and IDEAL. Breast density was calculated with semi-automated 
software. Statistical analysis was performed with non-parametric tests. 

Results: Mean percentage of breast density calculated in each sequence was: TI-tSE =56%; T2-tSE =52%; IDEAL FatOnly 
= 55%; IDEAL WaterOnly =53%, VIBRANT =55%. Significant differences were observed between T2-tSE and both T1-tSE 
(p<0.001), VIBRANT sequences (p = 0.009), TI-tSE and both IDEAL WaterOnly (p = 0.007) and IDEAL FatOnly (p = 0.047). 
Breast density percentage showed a positive linear correlation among different sequences: r>0.93. 

Conclusions: Differences exist between MRI sequences used to assess breast density percentage. Tl-weighted sequences 
values were similar to IDEAL sequences. 
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Introduction 

Breast density is a factor that should be taken into consideration 
when evaluating breast cancer risk. Breast density is an 
independent risk factor for cancer [1]: it has been demonstrated 
that women with dense breasts have a 4-6 times higher risk than 
women with fatty breasts [2—4]. A quantitative evaluation of breast 
glandular tissue is also very important for epidemiological studies 
[5,6]. Cancer risk predictive models, such as the Gail model, 
improve their predictive accuracy when breast density is incorpo- 
rated [7,8]. To date, there are different density classification 
methods like Tabar's classification [9,10], Wolfe's parenchymal 
patterns [11] and several quantitative evaluations of percentage 
mammographic density using automated or semi-automated 
computer-aided techniques [12-16]. However, the practical and 
logistical difficulties in correlating histopathological samples of the 
breast with radiological images prevent the definition of a gold 
standard for breast density evaluation. Therefore, a standard 
classification method to evaluate breast density is a critical 
challenge, in particular for Magnetic Resonance Imaging 
(MRI) [17, 18]. MRI has been used for breast density assessment 



[19] and correlation has been demonstrated between breast 
density percentage on MRI and mammography [17,20-24]. 
Other techniques, like mammography or tomosynthesis, could 
be adjusted on the basis of MRI for the evaluation of breast density 
percentage [25]. MRI assessment of breast density has been 
estimated using different sequences. Since 2007, sequences that 
can clearly separate the fatty non-glandular tissue from the true 
glandular tissue and its water content have been developed and 
made available for clinical purposes [26]. These sequences are 
called IDEAL (Iterative Decomposition of water and fat with Echo 
Asymmetry and Least squares estimation) and are believed to 
assess directly the biochemical features and composition of breast 
tissue similarly to physiology although they are still a surrogate of 
the real histological reference standard [26]. Other MRI 
sequences employed in the past evaluated breast density with 
proton density principles, as for mammography and digital breast 
tomosynthesis [25,27]. To our knowledge, there are no data in the 
literature comparing standard MRI sequences and IDEAL 
sequences for breast density percentage evaluation on a 3T MRI 
system. Therefore, the purpose of this study was to compare breast 
density percentage assessment with different MRI sequences 
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considering IDEAL sequences (Iterative Decomposition of water 
and fat with Echo Asymmetry and Least squares estimation) as a 
possible physiology-like reference. 

Materials and Methods 

Ethics Statement 

The study was approved by the Local Ethic Committee 
(National Institute for Cancer Research) and written informed 
consent was obtained from all participating women. 

Patients 

48 consecutive patients (mean age 41, years; range, 35-67 years; 
mean weight (kg) 56 ±4; mean BMI 22 ±3) underwent MRI 
examination, from March 2010 to October 2012. All women, if 
fertile, were examined between days 5 and 15 of their menstrual 
cycle. MRI examinations were performed after clinical and 
radiological examination, following EUSOMA (European Society 
of Breast Cancer Specialists) Guidelines [28] which were 
considered the inclusion criteria. Exclusion criteria were: breast 
prosthesis, bilateral breast pathology, claustrophobia, inability to 
tolerate MRI and lack of written informed consent. 

MRI 

MRI examinations were performed on a GE Signa HDx 3.0T 
scanner (General Electric Medical Systems, Milwaukee, WI, USA) 
with a dedicated eight-channel bilateral breast coil: the patient was 
placed in a prone position, without any compression of breasts. 
MRI scan protocol included the following standard sequences: Tl- 
tSE, T2-tSE, pre- and post-contrast agent (Gadobenate dimeglu- 
mine, MultiHance 0,5 M, Bracco, Italy) VIBRANT (Volume 
Imaging for Breast Assessment) with the adjunct of IDEAL 
(Iterative Decomposition of water and fat with Echo Asymmetry 
and Least squares estimation) sequences [25,29]. IDEAL sequenc- 
es are based on Dixon's method for fat suppression with the "in- 
and out-of-phase" technique: two images are acquired with 
different echo times (TE). Differences in chemical shift between 
fat and water should permit a clear separation of these tissues in 
perfect conditions [30]. This method was progressively modified 
and improved with the reduction of B0 field inhomogeneities 
[26,31,32]. IDEAL sequences were included in this protocol to 
distinguish fat tissue from water (FatOnly sequences and Water- 
Only sequences). WaterOnly sequences were considered to 
represent glandular tissue since fat is a water poor tissue. These 
sequences correlate well with physiology, but they are considered a 
surrogate for the real histological reference standard [26]. Main 
sequences parameters are reported in Table 1. 

Breast density analysis with semi-automated software 

Breast density analysis was performed only on the unaffected 
breasts. The data set comprised MRI images obtained in 48 
patients. Two patients were completely excluded for bilateral 
pathological findings. In the other 46 patients, the breast affected 
by cancer, fibroadenomas and cysts was excluded to avoid 
interference with the breast density analysis. There were 12 
patients screened for BRCA-1 mutations and had negative results 
after MRI, 25 patients with unilateral breast cancer (22 invasive 
ductal carcinomas and 3 invasive lobular carcinomas), and 9 with 
unilateral pre-malignant lesions. Every MRI sequence was 
evaluated before the administration of contrast agent to avoid 
possible influences of Gadolinium on breast density assessment. 
Breast density was calculated for each sequence using semi- 
automated software (MedDensity®). This software is a home- 
grown and previously validated for mammography, tomosynthesis 



and MRI [15,16,25]. All slices of every MRI sequence were 
included. We used the semi-automated version of the software 
(Fig. 1) because it permits proper tuning by the reader to guarantee 
more precise results. The software was adjusted to reduce user- 
related variability in assessing breast density with different 
sequences. The evaluations were performed by two radiologists 
independendy (M.C. and A.T.). The two radiologists had 22 and 7 
years respectively of experience in breast imaging and with more 
than 80 breast MRI reported every year. After the first evaluation, 
measurements were repeated after 2 months to assess intra- and 
inter-observer variability. 

Statistical analysis 

Statistical analysis was performed with commercially available 
softwares IBM SPSS Statistics v. 19 (International Business 
Machines Corp., New York, NY, US) and R (http://www.r- 
project.org/). After preliminary data checks, by mean of 
histograms, skewness and kurtosis index, the normality of 
distribution of each sequence, the non-parametric Friedman test 
was used to compare the five MRI sequences. If an overall 
significant result was detected, paired comparisons between the 
sequences were performed using the non-parametric Mann- 
Whitney test for paired samples. A p-value <0.05 was considered 
statistically significant. For multiple comparisons p-value was 
adjusted using the Holm correction criteria. Intra- and inter- 
observer variability were calculated using K statistics. Agreement 
was defined on the basis of Fleiss classification as follows: <0.40, 
poor; 0.40-0.59, moderate; 0.60-0.75, good; >0.75, excellent 
[28,33,34]. In addition, we considered in the analysis the effects of 
user interpretation and manual selection of threshold as 
confounding factors and adjusted the software interface and the 
statistical significance accordingly: p values <0.01 were considered 
statistically significant for that specific purpose. 

Results 

Forty-six women were included in the study. Mean percentage 
of breast density calculated in each sequence was: Tl-tSE = 56%; 
T2-tSE =52%; IDEAL FatOnly =55%; IDEAL WaterOnly 
= 53%, VIBRANT =55% (Table 2). Breast density percentage 
showed positive linear correlation among different sequences: 
r=0.95 between Tl-tSE and T2-tSE; r=0.93 between Tl-tSE 
and FatOnly; r=0.93 between T2-tSE and FatOnly (we reported 
the maximum and minimum values of r). 

Comparing the five MRI sequences, statistically significant 
differences were observed between T2-tSE and both Tl-tSE (p< 
0.001) and VIBRANT sequences (p = 0.009 using Student's T-test; 
p = 0.002 using non-parametric test). Complete results are shown 
in Table 3. 

Further statistically significant differences were observed 
between Tl-tSE and both WaterOnly (p = 0.007) and FatOnly 
(p = 0.047). Complete results are shown in Table 4. 

Intra-observer and inter-observer agreement of the two 
radiologists in the evaluation of breast density were considered 
to be very good (reader 1: k = 0.93; reader 2: k = 0.95; reader 1 vs 
reader 2: k = 0.95). 

Discussion 

Given the increasing importance of breast density percentage 
assessment in both the clinical and research fields, the develop- 
ment of a reliable and reproducible evaluation method would be 
useful [15]. Breast density percentage assessment has a very 
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Table 1. Sequence parameter settings. 



Sequence 


TR (ms) 


TE (ms) 


Flip angle 


Matrix 


Slice thickness (mm) 


Spacing (mm) 


Approximate Scan 
time (s) 


T2 


5200 


103 


90° 


350x350 


4.0 


3.0 


200 


T1 


600 


9 


90° 


350x350 


4.0 


3.0 


300 


VIBRANT 


6.2 


3 


10° 


350x350 


1.2 


1.2 


90 


IDEAL 


4380 


130 


90° 


360x360 


1.2 


3.0 


3600 



IDEAL :lterative Decomposition of water and fat with Echo Asymmetry and Least squares estimation. VIBRANT: Volume Imaging for Breast Assessment. 
s = seconds, mm = millimeters, ms = milliseconds 
doi:1 0.1 371 /journal.pone.0099027.t001 



relevant role in life-time cancer risk calculation [1—4], therefore its 
inclusion in predictive models may improve their accuracy [7,8]. 

To date, many methods for breast density assessment are based 
on the evaluation of analogic or digital mammography, but none 
of them are calibrated with a reference standard similar to 
histology. Using a method that is able to assess breast density 
percentage close to the real histological composition of the breast 



may be important to achieve a breast density percentage 
assessment similar to the real anatomy of breast tissue. Indeed, 
histology would give a true representation of the breast density, 
but histological specimens are generally not available for practical 
reasons. Recendy, MRI proved to be a potentially reliable 
technique for this purpose [27]. Three-dimensional MRI imaging 
of the breasts may be better than two-dimensional techniques 
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Figure 1. Quantitative assessment of breast density using the thresholding method. In this figure an example of the graphical computer 
interface is demonstrated. On the left, the edges of the breast are identified (edges are artificially thickened for visual purposes), then the radiologist 
adjusts the density threshold (semi-automated method) on the right, and finally the percentage of breast density is shown (black circle). In this 
example, the breast result was 41% dense. 
doi:1 0.1 371 /journal.pone.0099027.g001 
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Table 2. Mean, median, maximum and minimum values of breast density percentage evaluated on every MRI sequence with the 
semi-automated software. 





Sequence 


Mean 


Median 


Maximum 


Minimum 


Standard Deviation 


T2 


52 


58 


83 


6 


22 


T1 


56 


61 


95 


7 


23 


VIBRANT 


55 


56 


89 


8 


23.2 


WaterOnly 


53 


534 


95 


3 


23.1 


FatOnly 


55 


54 


92 


3 


23 



Results include both readings from both radiologists. IDEAL sequences (Iterative Decomposition of water and fat with Echo Asymmetry and Least squares estimation) 
are reported are WaterOnly and FatOnly. VIBRANT: Volume Imaging for Breast Assessment. 
doi:1 0.1 371 /journal.pone.0099027.t002 



similarly to other 3D imaging techniques that may be superior to 
2D evaluations [35]. In the literature, other authors evaluated 
breast density with MRI. Most of them used standard Tl- 
weighted sequences, yet it has not been proven whether or not this 
evaluation resembled the real composition of breast tissue. In 
particular, various researchers have used the following sequences 
for breast density assessment: Khazen et al. used pre-contrast Tl- 
weighted [17], Lee et al. used Tl -weighted spoiled gradient echo 
fast low-angle shot [22], Klifa et al. used 3D fat suppressed spoiled 
gradient echo pulse [19,36], Wei et al. used coronal 3D spoiled 
gradient recalled echo pre-contrast Tl -weighted [23], Thompson 
et al. used pre-contrast Tl -weighted [37], Eng-Wong et al. used 
Tl -weighted spoiled gradient-echo with fat suppression [38], and 
Wang et al. used a T 1 -weighted non-contrast fat-saturated images 
[39] . A common limitation of all these studies was that none of 
them used a sequence like the IDEAL that may be more likely to 
resemble the true breast tissue composition. These IDEAL 
sequences may be used as a physiology-like reference. A 
comparison between results obtained from MRI and mammo- 
graphy examination demonstrated that MRI and mammography 
had good correlation only in low-density breasts (BI-RADS 1 and 
2). In women with a greater breast density the correlation between 
MRI and mammography was low [36]. Evaluation of breast 



density on mammograms has limitations and technical problems. 
For this reason a three-dimensional imaging technique may be 
more accurate [40]. In the present study we added IDEAL 
sequences to study breast tissues from a biochemical perspective 
and not only according to the proton density typical of Tl weight 
[25,27]. This sequence may represent a reference in the 
assessment of breast density percentage because it is able to 
separate FAT and WATER. The main disadvantage of IDEAL 
sequence is longer scan time compared to other fat-suppression 
techniques [26]. Comparing IDEAL sequences with other MRI 
sequences in our study protocol (Tl- and T2-weighted and 
VIBRANT), we found that there were differences in breast density 
percentage assessment. Breast density percentage evaluated on 
T 1 -tSE and IDEAL WaterOnly (our reference standard for breast 
density) were similar. On the contrary, percentages evaluated on 
T2-tSE and VIBRANT sequences were different from IDEAL 
WaterOnly. T 1 -tSE showed similarity in the assessment of breast 
density percentage to IDEAL WaterOnly. From these data we 
suggest that Tl -weighted sequences produced images that are 
similar to the composition of breast tissues. The reliability of Tl- 
weighted sequences in the assessment of breast density is enhanced 
by our data. A Tl -weighted sequence is often included in the 
standard breast MRI protocol, while IDEAL sequences are not so 



Table 3. Test for paired samples to evaluate mean value differences between results of breast density percentage obtained on 
each MRI sequence. 





Sequence Comparison 


Differences in pairs 






P-value 




Mean 


SEM 


95% CI 












Inferior 


Superior 




T2 - Tl 


-4.25 


.98 


-6.22 


-2.27 


.000 


T2 - VIBRANT 


-2.85 


1.04 


-4.94 


-.75 


.009 


T2 - WaterOnly 


-1.46 


1.18 


-3.84 


.92 


.225 


T2 - FatOnly 


-2.00 


1.52 


-5.06 


1.06 


.195 


Tl - VIBRANT 


1.39 


1.24 


-1.10 


3.90 


.268 


T1 - WaterOnly 


2.78 


.90 


.97 


4.60 


.003 


T1 - FatOnly 


2.31 


1.19 


-.08 


4.71 


.058 


VIBRANT - WaterOnly 


1.39 


1.05 


-.73 


3.52 


.195 


VIBRANT - FatOnly 


.87 


1.33 


-1.82 


3.56 


.518 


WaterOnly - FatOnly 


-.55 


.60 


-1.76 


.66 


.365 



95% C.I. and P-value have also been calculated. IDEAL sequences (Iterative Decomposition of water and fat with Echo Asymmetry and Least squares estimation) are 
reported as WaterOnly and FatOnly. VIBRANT: Volume Imaging for Breast Assessment. 
SEM = Mean Standard Error, 95% CI = Confidence Interval for the difference. 
doi:1 0.1 371 /journal.pone.0099027.t003 
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widely available and increase the duration of the examination. For 
this reason we suggest that Tl -weighted sequences are sufficient to 
assess breast density on MRI. 

This study has several limitations. All examinations were 
performed on a GE Signa HDx 3.0T scanner but we suppose 
that for breast density assessment similar results are likely to be 
obtained on a lower magnetic field scanner. This assumption is 
supported by the fact that Tl -weighted sequences, for breast 
density assessment, are relatively stable at 1.5T and 3.0T. We 
acknowledge that further research could be performed to 
investigate this hypothesis. IDEAL sequences work better on a 
3.0T than on a 1.5T magnetic field because the chemical shift 
between water and fat is increased and consecutive echo groups 
have smaller spacing between them [26]. Image quality on 3.0T 
with a well-defined contrast among different tissues is improved. In 
our study we used semi-automated software because superiority 
control and regulations by reader was considered essential due to 
artifacts related to BO and Bl inhomogeneity [41]. A fully- 
automated software could have reduced time, costs and might 
have guaranteed more reproducible results, however an accurate 
control by the reader was sufficient to achieve similar results 
between different readings. Intra- and inter-observer agreement 
was very good and it is unlikely that the human reader may have 
influenced the results significantly. The software which was used 
was straightforward and easy to use and very little training was 
needed. No statistically significant differences between the first and 
the second reading performed after two months was found. For 
this reason radiologists' expertise in breast density assessment had 
litde impact. Concerning inclusion criteria we used EUSOMA 
Guidelines [28], therefore our patients represent a selected 
population of women for whom MRI was clinically suggested 
based on these guidelines. Institutions with different inclusion 
criteria may have a different patient population. This is the first 
study of this kind and our data were statistically significant even 
with a relatively small number of patients. Another limitation 
could be the artifacts related to B 1 inhomogeneity which should be 
considered when using a 3T MRI system. All the sequences were 
acquired consecutively on the same patient to reduce the 
measurement error due to Bl artifacts. Bl field homogeneity does 
not depend on body type and it is significantly improved by 
performing local RF shimming. Our MR system performs 
routinely local RF shimming for breast MRI. In the future, the 
value of a subject specific dual-transmit approach for improving 
Bl field homogeneity should be assessed for breast MRI. It has 
been demonstrated that with this approach, B 1 field homogeneity 
is significantly improved by performing local radiofrequency 
shimming with 2 independent radiofrequency-transmit channels 
[42] . In our study there was not a shimming with a dual-transmit 
system to improved local radiofrequency homogeneity, but only a 
standard shimming. Further research is needed to assess if a 
shimming with a dual-transmit system is able to influence breast 
density assessment on MRI. 

In conclusion we found significant differences among MRI 
sequences for breast density percentage assessment. If the IDEAL 
WaterOnly sequence is applied as a potential reference, the Tl- 
weighted sequences showed similar breast density percentage 
values as the reference. Therefore, in clinical practice, the Tl- 
weighted sequence could be a reliable and efficient sequence to 
measure breast density, and might possibly reflect the true 
composition of breast tissue. 
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